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We study the emergence of discontinuous shear-thickening (DST) in cornstarch, by combining 
macroscopic rheometry with local Magnetic Resonance Imaging (MRI) measurements. We bring 
evidence that macroscopic DST is observed only when the flow separates into a low-density flowing 
and a high-density jammed region. In the shear-thickened steady state, the local rheology in the 
flowing region, is not DST but, strikingly, is often shear-thinning. Our data thus show that the stress 
jump measured during DST, in cornstach, does not capture a secondary, high-viscosity branch of the 
local steady rheology, but results from the existence of a shear jamming limit at volume fractions 
quite significantly below random close packing. 


Granular suspensions are seemingly simple systems 
composed of non-colloidal rigid particles immersed in 
Newtonian fluids. Although the only relevant interac¬ 
tions are hydrodynamic and contact forces, they present 
a rich rheology [ 1 ], which includes shear-thinning or shear 
thickening, normal stress differences, shear banding and 
yield stress behaviour. With increasing packing fractions, 
continuous shear thickening (GST) appears at lower and 
lower strain rates and becomes more and more abrupt. 
This eventually leads to the emergence of a spectacu¬ 
lar phenomenon: discontinuous shear thickening (DST), 
an order-of-magnitude jump of the stationary (macro¬ 
scopic) stress when strain rate crosses some threshold 
value. DST is observed in a narrow range of packing frac¬ 
tions near random close packing and recent works suggest 
that it is a dynamic or shear jamming transition. This 
motivates a considerable attention from experiments [ 2 - 
4], numerical simulation [5-7] and theory [ 8 , 9]. 

The origin of shear thickening is essentially unknown. 
Proposed mechanisms include: the crossover to Bagnold 
scaling due to particle inertia [ 10 ]; the coupling between 
normal and shear stresses associated with material dila- 
tancy [11, 12]. In numerical simulations DST was ob¬ 
served in systems which are forced to remain homoge¬ 
neous by periodic boundary conditions [5]. It was thus 
suggested that it might arise from a jump of the local 
response, which occurs at fixed packing fraction, and is 
caused by the proliferation of frictional contacts [5, 9]. 

Examples abound, however, in granular suspension 
rheology where the macroscopic response is sharply dif¬ 
ferent from the local behaviour: the apparent yield 
stress [13, 14] or a transient DST behavior [10], for ex¬ 
ample, were both shown to result from the emergence 
of flow inhomogeneities due respectively to density dif¬ 
ferences and migration. The observation of finite size 
effects in DST measurements [11, 12, 15], at sizes up to 
100 particles, thus questions the idea that DST is a purely 
local phenomenon. In this context, it is unclear how a 
phase diagram for the local response can be deduced from 
macroscopic experimental data and compared with nu¬ 
merical [5], and theoretical [9] works. 

This situation motivates us to study DST in corn¬ 


starch [16-18], the paradigmatic system for this phe¬ 
nomenon, by accessing experimentally the local response. 
To this aim, we have designed a new velocity-controlled 
rheometer that can exert steady torques up to 10 N.m 
on a wide-gap Couette cell inserted in our Magnetic Res¬ 
onance Imaging (MRI) scanner. It enables us to access 
local velocity and particle volume fraction profiles in the 
flowing sample [13, 19], for a wide range of strain rates 
( 7 ) covering GST and DST. We thus find that macro¬ 
scopic DST is observed when the flow localizes, i.e. the 
system separates into a low-density flowing and a high- 
density jammed region. In the shear-thickened steady 
state, the local rheology in the flowing region is not DST 
but, strikingly, is most often shear-tbinning. Our data 
support that DST should be attributed to the existence 
of a shear jamming limit 7 dst( 0 ) at volume fractions (j) 
significantly smaller than random close packing. 

The cornstarch (Sigma) is the same as in [ 11 , 16]: it 
presents slightly polydisperse particles of size ^ 20 /im 
with irregular shapes. We have measured the random 
close packing fraction 0 rcp — 55% by tapping the pow¬ 
der for long time in a container [20]. We have also mea¬ 
sured the random loose packing (/)rlp — 35% as the 
density achieved after sedimentation of a dilute suspen¬ 
sion [21, 22]. The suspensions are prepared by mixing 
cornstarch in demineralized water at initial volume frac¬ 
tions ranging from 33.5% to 43.9%. 

Let us start with macroscopic rheometry measure¬ 
ments using a Malvern Kinexus Pro stress-controlled 
Couette rheometer. The inner and outer cylinders radii 
are: Ri = 12.5mm and Rq = 18mm; the inner cylin¬ 
der height is h = 37.5mm; the gap is large enough to 
avoid confinement effects [11, 12, 15]. Both cylinders are 
roughened to avoid wall slip. We report on Fig. l-(a) the 
torque T measured on the inner cylinder during a 3 min 
logarithmic ramp of Ft (the inner cylinder rotation rate) 
from 5.10“^ rpm to 2.10^ rpm. Near the inner cylinder, 
the local stress is 

T{Ri,n)=T{n)/{27rR^h) ( 1 ) 

the local strain rate is a priori unknown, but is classicaly 
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estimated as: 


^{Ri,n)c^2nRl/{Rl-Rj) ( 2 ) 

Combining Eq. ( 1 ) and ( 2 ) yields the apparent viscosity 
77 = r/y vs 7 relation reported on Fig. l-(b). Note that 
at low T] macroscopic inertial effects may arise (Taylor- 
Couette instabilities), which limit the accessible 7 range 
at low (/)’s. 

At low 7 , the response is clearly shear-thinning. It 
crosses over to CST around a characteristic strain rate 
7csT defined at that where 77 reaches a minimum. Vis¬ 
cosity jumps at some higher strain rate 7dst which char¬ 
acterizes the onset of DST. 

The classical estimate of ^{Ri) as Eq. (2), which is 
based on the Newtonian solution, remains unsatisfactory. 
An exact expression exists so long as the material is ho¬ 
mogeneous in the studied torque range: 


W;To) = 2^ (t— 


n=0 


dCl 




(3) 


It cannot be used around the DST transition, however, 
because the apparent singularity of dT/dO cannot be re¬ 
solved experimentally. Nevertheless we have checked that 
using the first two terms of (3) gives the same qualitative 
behaviour as in Fig. l-(b) up to DST. We will show be¬ 
low that the flow does remain homogeneous below DST, 
which allows us to use (3) to accurately compute the 
7csT and 7dst values that are reported on Fig. l-(c). 
Interestingly, (i) at any volume fraction CST is always 
observed before DST; (ii) both ycsT and 7dst seem to 
vanish around the same volume fraction 0 c which we es¬ 
timate to be ^ 45% <C 0RCP by linear extrapolation. 

Now we turn to velocity-controlled MRI rheometry. 
Our Couette cell has inner and outer radii Ri = 3cm and 
Ro = 5cm (resp.) and inner cylinder height h = 11cm. 
Both cylinders are roughened to avoid slip, which we 
checked from velocity profiles. All experiments discussed 
below are performed by preparing a homogeneous mate¬ 
rial with mean volume fraction 0o = 43.9% (experiments 
at 40%, 41% and 42.5% show similar features). Note 
that our MRI cell dimensions differ from those of the cell 
used to obtain the macroscopic data of Fig 1 . Hence, to 
capture DST with this new cell, the range of rotation ve¬ 
locities must be adapted: at 0o = 43.9%, we use ft values 
between 5rpm and 100 rpm. In each MRI experiment, O 
is held fixed until steady state is reached. 

Our MRI [13, 23] provides the stationary packing frac¬ 
tion 0 (r, Q) and azimuthal velocity v{r, D) at any radial 
position r. From the latter we extract the local strain 
rate j(r,Q) = vjr — dvjdr. We do not have access 
to torque measurement. But, since the local stress is 
r(r, Vt) = r{Ri,fl)R‘f/r‘^ in the Couette geometry, we can 
estimate the local viscosity profile as: 


ri{r,n) 


T{Ri,n) Rj 
7 (r, r 2 


(4) 
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FIG. 1. (Color online). Macroscopic rheometry data, (a) 
Torque T vs rotation rate Q during a logarithmic ramp at 
various packing fractions 0. (b) Apparent viscosity vs ap¬ 

parent shear rate extracted from (a), (c) Critical shear rates 
7CST and 7 dst vs 0. (d) Reversibility test at 0 = 43.9%: 
succession of up, down, and up ramps. 


up to the unknow prefactor r(Ri, f 2 ). 

In Fig 2, we plot the steady velocity and concentration 
profiles thus measured for a few f2s. Velocity is normal¬ 
ized by its value at the inner cylinder. 

For small f^’s (5 and 7 rpm), we find that the density 
profile is homogeneous while the flow extends through¬ 
out the gap. It is known that non-Brownian suspensions 
may slowly become inhomogeneous due to shear-induced 
migration and sedimentation [13, 14]. We checked that 
this does not occur before strains larger than a few thou¬ 
sands, which is much larger than the strain range over 
which we collect data. Since density is uniform, we can 
access the local rheology 77 ( 7 , 0o) as follows. For each 
matching the single parameter T{Ri^Q) in Eq. (4) pro¬ 
vides the complete 77 (r, ^4) profile, which can thus be plot¬ 
ted vs 7 (r, f2). These values are compared on Fig. 3-(a) 
with the local rheology data obtained from macroscopic 
rheometry using Eqs. (1) and (3) to evaluate the local 
stress and strain rate near the inner cylinder. Clearly, in 
the shear-thinning regime, the stationary local response 
within the gap matches our macroscopic measurements 

This local homogeneous response is observed so long 
as the maximal local strain rate, which is reached at the 
inner cylinder, lies below the 7dst(0o) value identified 
in macroscopic rheometry [Fig. l-(c)]. Homogeneity and 
locality then enable us to estimate the critical rotation 
rate at which DST is expected as: Odst(0o) — 8 rpm. 

A sudden transition occurs as soon as crosses 
^dst(0o)- As shown on Fig. 2-(a), the flow then abruptly 
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FIG. 2. Steady MRI data for a 0o = 43.9% cornstarch suspen¬ 
sion and different rotational velocities Q. (a) velocity profiles, 
(b) density profiles; solid lines indicate 0rlp and 0rcp- 


stops in a large region. Namely, the velocity profile jumps 
from one of the rightmost curves, corresponding to homo¬ 
geneous flows, to the leftmost one, i.e. the most strongly 
localized flow. Note that measuring a single velocity pro¬ 
file requires accumulating MRI data over ^ 30s, which 
corresponds here to a strain of order 50. Upon cross¬ 
ing Udst(0o): the first measurable velocity profile is al¬ 
ready localized. The flow subsequently remains steady 
over thousands units of strains. DST is thus clearly con¬ 
comitant with shear localization. 

As Q increases further, the velocity profiles progres¬ 
sively extend to the right (i.e. towards the outer cylin¬ 
der). In all cases, the system remains separated into a 
flowing layer near the inner cylinder and a jammed re¬ 
gion near the outer one. The fraction of the gap that is 
jammed is reported on Fig. 3-(b): it jumps at Udst((?^o) 
and then slowly decreases. 

Comparing these velocity profiles [Fig. 2-(a)] with lo¬ 
cal density data [Fig. 2-(b)], we find that, quite strik¬ 
ingly, the flow localization at U = lOrpm ^ Udst((?^o) 
is associated with the sudden emergence of density in- 



FIG. 3. (a) Gomparison of local rheometry data obtained 

from MRI measurements (open symbols) and near the inner 
cylinder in macroscopic rheometry (filled symbols) in homo¬ 
geneous conditions. Upper data: 0 = 43.9%; lower data: 
0 ~ 33.5%. (b) From the velocity profiles [Fig. 2-(a)]: frac¬ 
tion of the gap which is jammed vs U. 


homogeneities. Namely, the volume fraction decreases in 
the flowing layer, while it increases in the jammed region, 
as required by the conservation of particle number. As 
Q increases beyond Udst(0o): the progressive extension 
of the flowing layer is accompanied by a broadening of 
the low-density region. At high strain rates, the den¬ 
sity saturates, in the flowing layer, at a packing fraction 
0min — 33% ^ 35% 0RLP and, in the jammed region, 
at 0 ^ (j)RCP • It is noteworthy that the density profile can 
achieve multiple forms depending on shear history [24]. 

Let us emphasize that the change of density created by 
the DST transition is irreversible. Indeed, once a station¬ 
ary profile 0(r, Ui) has been produced by ramping U up 
to some arbitrary > Urst, we found that the density 
profile remains the same under any subsequent lowering 
of U. This irreversibility shows up in our macroscopic 
rheometry setup (the small Couette cell) as illustrated 
on Fig. l-(d) where we plot the torque T vs U during: 
(i) an initial up-ramp that drives the system through 
the DST transition; followed by (ii) a down-ramp. The 
torques measured during the up- and down-ramps clearly 
lie on different branches. However, if we subsequently 
(iii) re-increase U, torque T(U) tracks the data previ¬ 
ously obtained on the down-ramp. Hence we reason that 
on the down (ii) and up (iii) ramps, the system explores 
quasi-reversibly a family of states associated with some 
inhomogeneous density profile set by Ui, the maximum 
value of the rotation rate achieved during shear history. 

We now probe the local rheology within the shear- 
thickened flow. Three of our datasets, for U = 40, 60, 
and lOOrpm, present a finite-width region of roughly uni¬ 
form packing fraction (j)min — 33%. Using only the data 
collected in these regions, we can hence implement the 
same approach as in the fully homogeneous case, to ac¬ 
cess the local rheology r]{j,(j)) at (p = 0min- For each U, 
we match the prefactor r(i?i, U) and plot in Fig. 3-(a) the 
resulting rj{r, Q) vs 7 (r, U) data points. Clearly, the local 
0min) rheology is shear-tbinning, and extrapolates 
the 0 = 33.5% rheometry measurements. 
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FIG. 4. Phase diagram: scatter plot of the local 
( 7 (r, Q), 0(r, Q)) values compared with the 7 cst( 0) (red) and 
7dst(0) (blue) lines obtained from macroscopic rheometry 
[Fig. l-(c)]. 


These data demonstrate that macroscopic DST cannot 
be straightforwardly interpreted as capturing a compara¬ 
ble change of the local response at fixed volume fraction. 
Instead, it results from the separation of the system into a 
low-density-flowing and a high-density-jammed regions. 
Somewhat paradoxically, the local response was found 
above to be shear-thinning in the flowing layer of the 
macroscopic shear-thickened state. 

In the flowing layer, the density plateau only devel¬ 
ops around 0min- It is hence at this packing fraction 
only that we can extract the local rheology 77 ( 7 , 0) from 
MRI data. To further qualify the local response, we 
now construct a phase diagram for all the local flow and 
jammed states. To this end, we collect the operating 
points ( 7 (r, f^), 0(r, f7)) obtained with the MRI for all f^’s 
and all r’s and report them on Fig. 4, where they are com¬ 
pared with the curves 7 cst( 0) and 7 dst( 0) that delimit 
the flow regimes previously identified using macroscopic 
rheometry. Our MRI data include points obtained in the 
flowing layer as well as in jammed regions: in the scatter 
plot, they are respectively found on the left (low 0, with 
non-vanishing 7) and right (high 0, 7 ^ 0) hand sides. 

Local flow point fall in either the shear-thinning or 
CST domains, which supports that both types of lo¬ 
cal rheologies are found. Moreover, they all lie below 
the 7 dst( 0) line, which thus appears to be an intrinsic 
material limit beyond which no steady flow can be sus¬ 
tained. No flow is found in particular at packing fractions 
above the value 0c — 45% where 7cst(0) and 7dst(0) 


vanish. As 0c is quite distinct from the jamming point 
0RCP — 55%, a continuum of jammed states is accessible 
beween 0c and 0rcp- 

In the literature, experimental observations of macro¬ 
scopic DST are usually analyzed while assuming that the 
flow remains homogeneous: the DST stress jump hence 
is interpreted as picking up a secondary, high-viscosity, 
branch of the steady local response. As we have shown, 
however, homogeneity does not hold. The stress jump is 
a property of the global response, and is associated with a 
global reorganization of the flow, while the local response 
always lies within the shear-thinning or CST regimes. Of 
course, this is shown here for cornstarch, but we note that 
today, there is not in the literature, a single experimental 
data point that can be trusted to pick up the presumed 
secondary branch of s-shaped rheologies, since there is no 
other local observation of the flow response. 

Moreover, the DST scenario that we are unveiling here 
is rather unexpected in view of recent models based on 
s-shape flow curves [9]. Indeed, at the onset of macro¬ 
scopic DST, the part of the flow that reaches locally a 
limit of stability, i.e. the 7dst line, is the region near 
the inner cylinder, where the local shear rate is maximal. 
From s-shape models, in their current formulation, we 
would expect this locally unstable region to jump to a 
high-viscosity state. In stark contrast, we find that when 
DST occurs, this region jumps toward a low-density, low- 
viscosity shear-thinning state. Meanwhile jamming oc¬ 
curs in the region near the outer cylinder which, at the 
onset of DST, was well within its stability limit. 

It finally appears that DST is determined by two ma¬ 
jor factors: (i) the existence of a limit 7dst(0) beyond 
which no steady flow can be accessed; (ii) the fact that 
as 0c ^ 0 RCP 5 a continuum of jammed states exists. 
Indeed, due to (i) the inner cylinder region eventually 
reaches the 7dst(0) shear-jamming limit, that it can 
only escape by decreasing its local volume fraction; it 
thus drives migration throughout the gap, because (ii) 
the growing jammed layer can store material significantly 
above 0c. Both (i) and (ii) are direct consequences of a 
unique physical phenomenon, shear jamming. 

It is crucial here that the 7dst(0) line is not vertical, 
hence can be crossed by increasing 7 at given 0. It is 
so because the asymptotic flow regimes (high and low 7, 
resp.) have two different limit volume fractions (t)RLP 
and 0c, resp.). This contrasts sharply with [10] where 
jamming is rate-independent and DST purely dynami¬ 
cal. It remains open to identify the microscopic mecha¬ 
nisms that account for the splitting of the jamming limit 
observed here. Several may be envisionned, such as the 
increased mobilization of frictional contacts with increas¬ 
ing shear rate [5], or the angularity of particles, which is 
significant for cornstarch, or finally particle inertia, which 
might be causing the CST [10] that always preceeds DST. 
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